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Abstract 

We study the pinning quantum phase transition in a Tonks-Girardeau gas, both in equilibrium 
and out-of-equilibrium, using the ground state fidelity and the Loschmidt echo as diagnostic tools. 
The ground state fidelity (GSF) will have a dramatic decrease when the atomic density approaches 
the commensurate density of one particle per lattice well. This decrease is a signature of the 
pinning transition from the Tonks to the Mott insulating phase. We study the applicability of the 
fidelity for diagnosing the pinning transition in experimentally realistic scenarios. Our results are in 
excellent agreement with recent experimental work. In addition, we explore the out of equilibrium 
dynamics of the gas following a sudden quench with a lattice potential. We find all properties of 
the ground state fidelity are reflected in the Loschmidt echo dynamics i.e., in the non equilibrium 
dynamics of the Tonks-Girardeau gas initiated by a sudden quench of the lattice potential. 

PACS numbers: 03.75.Kk, 05.30.-d, 03.65.Yz, 67.85.De 
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I. INTRODUCTION 



In the past two decades, ultra-cold atomic systems have emerged as ideal playgrounds for 
the controlled simulation and manipulation of textbook models from many-body physics 
Using the full armory of developed tools the parameters of the underlying Hamiltonian can 
be tuned with an unprecedented precision allowing for the exploration of phase diagrams 
synonymous with condensed matter physics. In addition, the high degree of isolation, tun- 
ability and long coherent time scales associated with ensembles of ultra-cold atoms allow for 

n 

excellent time resolution of quantum dynamics |2J. 

For a long time studies of integrable systems were considered a purely academic pursuit, 
but by now can be created in the laboratory with ensembles of cold atoms. By applying 
the appropriate lasers to Bose-Einstein condensates, one dimensional arrays of atoms may 
be formed In the limit of strong interactions these arrays were observed to be in a 
fermionised state known as the Tonks- Girardeau gas 4|, |5j, a prototypical integrable system. 

In this work we will focus on the Tonks- Girardeau gas in a particularly interesting con- 
figuration which admits critical point. If a weak periodic potential is applied along the axial 



direction of a one-dimensional ultra-co 
simulation of the Sine-Gordon model 



d quantum gas it is possible to generate an atomic 
6J. When the interactions between the particles in 
the gas are sufficiently repulsive and the lattice is commensurate with the particle density 
(one particle per lattice well) this model has a quantum phase transition (at T OK) 
where atoms become 'pinned' to the Mott insulator state. In contrast to the well known 
superfluid-Mott insulator transition, pinning to the Mott phase occurs for an infinitesimally 
weak lattice potential js]. A spectacular recent experiment demonstrated this transition for 
ensembles of one dimensional ultra-cold gases 

In general, a quantum many-body system which undergoes a quantum phase transition 
may be written as 

H(X) = H + XH>, 

where A and H' are the driving parameter and the Hamiltonian driving the quantum phase 
transition (QPT) respectively. A feature of a phase transition is that if the parameter A 
is varied across the critical point, the energy spectrum undergoes a dramatic change i.e., 
the ground states of H(X) and H(X + 5X) will significantly differ. As a consequence the 
overlap of the ground states is expected to be sensitive to a QPT 8]. According to the 
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To nks- Girardeau gas has a pinning quantum phase transition at A = when the driving 
Hamiltonian (H 1 ) includes an optical lattice Vi(x) = Visin 2 (kx) commensurate with atomic 
density where the amplitude of the lattice is the parameter driving QPT i.e., V\ = 5X. In this 
paper we shall denote H(\ = 0) as the Hamiltonian of TG gas in a trapping potential Vq(x) 
and H(X + 5X = V{) as the Hamiltonian of TG gas in Vq(x) + Vi(x) potential. We denote 
the ground state of H{0) as \^>o) and ground state of H(Vi) as |$o)- We expect that the 
overlap of ground states (^ol^o) will be sensitive even to infinitesimally weak optical lattice 
if lattice periodicity is commensurate with atomic density {fil, 8(. In quantum information 
theory, the square modulus of the overlap is known as fidelity |9( and is a central concept in 
state characterisation. The ground state fidelity is defined as 

F=\(V \%)\ 2 . 

In this work we use this fidelity to study pinning quantum phase transition in the Tonks- 
Girardeau gas. We find, as expected, that GSF decreases with the increase of the lattice 
amplitude and size of the system. We emphasize that in the thermodynamic limit the GSF 
can unequivocally determine the pinning quantum phase transition for an infinitesimally 
small lattice amplitude. Nevertheless, the auxiliary trapping potential and finite size effects 
are important for experimentally relevant numbers of particles. We find that the GSF 
is in agreement with recent experiments on the pinning quantum phase transition (QPT) 



for a Luttinger liquid of strongly interacting bosons 



7|. All of the observed properties 



of ground state fidelity are also reflected in the dynamical evolution of the system i.e in 



survival probability or the Loschmidt echo 



lCH14j (for a review see e.g. 15|). The average 



value of the Loschmidt echo decreases for lower value of ground state fidelity; that is a 
general observation. Details of Loschmidt echo dynamics, such as the dominant frequency 
of revivals, depend on the particular trapping potential. We find that for the TG gas in 
an infinitely deep box potential oscillations of the Loschmidt echo are large and occur with 
smaller frequency in the critical region than in rest of parameter space. In the harmonic 
oscillator potential, the frequency of the Loschmidt echo revivals is constant until we reach 
a critical number of particles N P i nn , and after N P i nn the oscillations become irregular. 
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II. THE PINNING TRANSITION IN A TONKS-GIRARDEAU GAS 



Consider a gas of bosons confined in a tight waveguide at T « OK temperature with tight 
transverse trapping frequencies such that u± ^> fi/h, where \i is the chemical potential. 
In this regime we may describe the many-body system by an effective one dimensional 
Hamiltonian, 

H = [ dx&(x)[—4^ + V (x)]if(x) + - [ dx&(x)&(x)V(x)if(x), (1) 
J 2m ax 2 2 J 

where Vq(x) is an arbitrary one dimensional longitudinal external potential and g describes 
the strength of a short ranged interaction. In such one dimensional systems it is typical to 
introduce the following dimensionless parameter, 7 = mg/(h 2 p), which is the ratio of the 
kinetic energy to the interaction energy (p is the linear density). In the spatially uniform 
case the spectrum is gapless for all 7 and described by a Luttinger liquid of bosons. Let us 
assume a one dimensional optical lattice Vi(x) = Vi sin 2 (/cx) is applied along the longitudinal 
direction of the waveguide in addition to already existing trapping potential Vq(x), in this 
case Vi is the strength of the applied lattice and we introduce wave vector k = 2ir/\. When 
interactions are weak, 7 <C 1, and the lattice strength is much larger than the recoil energy 
Vi ^> Ejt = (hk) 2 /(2m), Eq. (CD) may be mapped on to the Bose-Hubbard model in the tight 
binding approximation [lj]. In this model, there is a phase transition as one changes the 
ratio of tunneling to atom-atom interactions, between a superfluid state where the atoms 
are free to tunnel between the wells coherently and a Mott state with an excitation gap and 
fixed number of particles per lattice site. 

Interestingly, in the opposite case when the strength of the applied lattice is much smaller 
than the recoil energy V\ <C Er, the Bose-Hubbard model is not applicable as the bosons 
now occupy several vibrational states in each well. In this case it was shown by Biichler 

n 

et al that the system maybe mapped to the famous Sine-Gordon model [6|, an effective 
low energy theory has been extensively studied in the literature as a rare example of an 
exactly solvable quantum field theory. Biichler et al showed that when the gas is in strongly 
interacting Tonks Girardeau limit, 7 ^> 1, and the lattice is commensurate with the density 
then the system will be 'pinned' to the Mott insulator state for an arbitrary weak lattice |6|. 
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III. THE FERMI-BOSE MAPPING THEOREM AND GROUND STATE FI- 
DELITY 



The pinning phase transition is quite straightforward to understand in the Tonks- 
Girardeau limit of strong repulsive interactions, g — > oo, in which this work will focus 
on. Physically, one may understand the pinning phase transition in this limit as the com- 
petition between the average inter-particle distance due to the strong interactions and the 
period of the potential. In this limit the hard core interactions play the role of the Pauli 



exclusion principle and the Fermi-Bose mapping theorem of Girardeau applies [16j. This 
theorem proves that the wavefunction of the system defined by a Hamiltonian such as Eq. ([T]) 
with g — > oo is equivalent to the properly symmetrised wavefunction of a gas of noninter- 
acting fermions in the same trapping potential Vq(x). As is customary for non- interacting 
fermions with periodic boundary conditions, an applied commensurate lattice Vi(x) leads to 
the opening of a single particle band gap of width A = Vj/4. This is the Mott insulating 
phase. 

As we will focus on the pinning transition in the Tonks Girardeau limit, let us briefly 
review the Fermi-Bose mapping theorem. The essential idea is that one can treat the inter- 
action term in Eq. ([I]) by replacing it with a boundary condition on the allowed manybody 
bosonic wave-function 

y B (xi,X2,...,x n ) = if \x i -x j \=0, (2) 

for % 7^ j and 1 < % < j < N. This is a hard core constraint meaning no probability exists 
for two particles ever to be at the same point in space. 

This constraint is automatically fulfilled by the corresponding noninteracting fermionic 
system using a Slater determinant such that 

1 N 

x 2 , . . ■ , x N ) = —=- det [ip n (xj)] , (3) 
ViV! n >i =1 



where the ip n (x) are the single particle eigenstates of the noninteracting system in trapping 
potential Vq(x). This, however, leads to a fermionic rather than bosonic symmetry, which 



can be corrected by a multiplication with the appropriate unit antisymmetric function 16| 



^b(xi,x 2 ,...,x n )= sgn(x i -x j )^ F (x 1 ,x 2 ,...,x N ), (4) 

l<i<j<N 



The power of the mapping theorem is that certain important many-body quantities of the 
To nks- Girardeau gas in an arbitrary external potential, can now be calculated using single 
particle states. The analytic nature of the many-body states of the gas in this limit are 
convenient to explore the properties of the pinning transition. 

A feature of the pinning quantum phase transition is that even a weak lattice can change 
the energy spectrum dramatically and the overlap of two ground states decreases. Using FB 
mapping, the ground state fidelity can be expressed via single particle basis 

N N 
01 i=l a 2 j=l 

1 N 

= iM^B-)"(-r 



(71 CT2 %=1 

|2 



i=l J 

detA| 2 (5) 



where elements of matrix A are 

A i:j = J '4)*(x)(f) j (x)dx. (6) 

If the system is in the ground state |\l/o) an d we suddenly turn on optical lattice V[(x), the 
probability that we will excite the s yste m away from the initial ground state is conveniently 



related to the ground state fidelity [18| 



p exc = i-K^oi$o>r. (7) 

In section [VI] we explore non-equilibrium dynamics after a sudden quench of lattice am- 
plitude. The fidelity of the TG gas is formally equivalent to a gas of non-interacting spin 



polarized fermions 191 ] . 



IV. PINNING TRANSITION FOR A TG GAS IN AN INFINITELY DEEP BOX: 
GROUND STATE FIDELITY 

In this section we apply the concept of ground state fidelity (GSF) [Eqs. ([5]) and ([6])] to 
study the pinning quantum phase transition for a TG gas in an infinitely deep box potential, 

fO, if < x < L 
(8) 
oo, otherwise . 
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The lattice potential is denned as Vi(x) = Vicos 2 (kx + 0). The periodicity of the lattice 
corresponds to the length of the box, k = Mir/L, where M is an integer. Thus, for = 
we have exactly M wells within the box, and we expect to see the signature of the pinning 
transition at N = M, where N is the number of particles. For the the other phases 
there are M — 1 well defined wells, and two half-wells at the edges of the box. In the 
thermodynamic limit the differences due to boundary effects will disappear (or become 
irrelevant), however in our simulations we will investigate these finite size effects which can 
be relevant for experiments. First we study the GSF numerically as a function of number 
of particles N for different lattice amplitudes V\ and different system sizes L. 

A. Numerical simulations 

In our numerical simulations the x-space grid is in units xq = \\ira. The lattice amplitude 
V\ and all other energies are in units of the recoil energy Er = (hk) 2 /(2m). The mass m 
corresponds to rubidium atoms 87 Rb. We shall fix the lattice wave vector to be k = 4tiXq 1 /3 
(A = 2n/k = 3xq/2), and keep it constant throughout this section. The length of the 
box L = Mir/k = MX/2 will vary. In all simulations Vi < Er, i.e., we are in the weak 
lattice regime Single particle (SP) states of Vq(x) are ip n (x) = \/2/ L sin(ri7rx / L) 

(n = 1,2,3,...). The SP states <f) n (x) of Vo(x) + Vi(x) are calculated numerically. From 
these one obtains GSF via Eq. flHJ). 

Figure Ufa) shows the GSF as a function of the number of particles N, for different values 
of the lattice amplitude. The phase = 0, i.e., Vftx) = V\ cos 2 (/cx). The size of the box is 
L = 40ir/k = 30xo, that is, M = 40 and we expect the pinning to occur at N = 40. Indeed 
we observe a dramatic decrease of fidelity when approaching commensurability, however, the 
GSF is equal for N = 39 and N = 40. One can argue that in the thermodynamic limit there 
is a single point at which the pinning takes place, and that this anomaly is a consequence of 
finite size effects. Nevertheless, such finite size effects are important for experimental systems 
as they occur at the relevant densities. The aforementioned anomaly will be explained in 
the next subsection using first order perturbation theory. We point out that the GSF 
obtained with the first order perturbation theory (black crosses) for N < M + 1, developed 
in Subsection IIV B\ is in perfect agreement with numerics for small amplitude V\ = O.IEr 
(blue circles in Fig. HJa)), for higher amplitudes there are discrepancies between first order 
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FIG. 1: Ground state fidelity (a) as a function of the number of particles N, for constant box 
size L = 30xo an d lattice wave vector k = Attx^ 1 /3 with different lattice amplitudes V\. The 
black crosses are analytical results (for N < M + 1) obtained with first order perturbation theory. 
Ground state fidelity (b) as a function of N, for constant lattice amplitude Vi = 0.3Er and wave 
vector k = Att/3 with different box sizes L. See text for details. 

perturbation theory and exact numerics in the dip of GSF, while outside of the GSF dip 
agreement is fairly good for all amplitudes, see Fig. HJa). 

Figure CO^b) shows the GSF as a function of the the number of particles N, for different 
values of L (lattice amplitude is held constant at the value V[ = 0.3Er). We clearly see that 
GSF decreases in the region of criticality with the increase of L, as expected. In the next 
subsection we will show that at the critical point F — > as L — > oo. 

Let us discuss the boundary effects for a finite-size system. Interestingly, if we use the 
phase (p = 7r/2 for the lattice, such that Vftx) = Vi sm 2 (kx), we obtain approximately 
identical values for the fidelity. In Fig. [2] we show GSF for cosine squared (0 = red circles) 
and sine squared (0 = n/2 blue crosses) lattice, these values overlap and come in pairs. This 
symmetry is lost for phases in between and ir/2. As an example, Fig. [2] shows GSF 
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as a function of the number of particles for = tt/3 (green squares) and = it/ 8 (pink 
triangles); there is a single point at which GSF has a minimum, either at N = 39 or at 
N = 40. For other phases, in between = and = tt/2, results are qualitatively and 
quantitatively similar as for = ir/3 and = tt/8. 




FIG. 2: Ground state fidelity as a function of N, for different phases of the lattice Vi(x). Outside 
presented interval of N, GSF is approximately identical for all phases. See text for details 



It is interesting that for a cosine squared lattice with exactly M wells the signature of the 
pinning occurs at N = M, M — 1, and that the same values are obtained for a sine squared 
lattice with M — 1 wells plus two half- wells. In order to see the differences between different 
boundary conditions we need to look at another quantity rather than the GSF. We choose to 
investigate the behavior of the energy (following the experiment 3]), and the single-particle 
density. 

In Fig. [3] we plot the single-particle (SP) energy spectrum of the potential potential 
Vq{x) + Vi(x) for both sine- and cosine- squared lattice. As in Fig. Ufa) the parameters are 
L = 30xo and k = Attxq 1 /3 (which gives M = 40), and the lattice amplitude is Vi = 0.85-Er. 
We see that the energy spectrum is different for these two lattices. Even though we cannot 
strictly speak about a gap for a finite size lattice, by observing Fig. [3] we see that the gap- 
like opening in the spectrum occurs at n = 40 (n is the index of a single-particle state) 
for a cosine-squared lattice, whereas for the sine-squared lattice it occurs at n — 39. These 
signatures for the pinning transition are intuitively expected when we think of the number of 
particles versus the number of wells in these two lattices. Even though we speak here about 



the SP spectrum, the energy gap wil 



be present in many-body excitations of the TG g clS clS 



well, because of the FB mapping 16|. We emphasize that 'gap' in SP spectrum occurs for 



9 




FIG. 3: Single particle SP energy spectrum for cosine-squared (f) = (black crosses) and sine- 
squared 4> = it/2 (red circles) lattice, other parameters are L = 30xo and Vi = 0.85Er. See text 
for details. 

SP states with the same (or approximately the same) wavelength as the lattice wavelength; 
this fact will be used in Section [V] 

We now turn to the single-particle densities p(N, x) = J2n=i \4>n(%)\ 2 / N , which are plotted 
in Fig. H] for the cosine- (a) and the sine-squared (b) lattice, in comparison with the density 
Po(N,x) = Yln=i li^nix)] 2 / N , and the lattice maxima and minima (parameters are N = 
M = 40, Vi = 0.55Er and L = 30xo). We see that in both cases the density maxima occur 
at the lattice minima as expected; the two cases differ at the boundary which is reflected in 
the energy spectrum but not in the GSF. 

In Fig. (a) we show the inset of the density p(N = M, x) (red crosses) vs. two nearby 
densities p(N = M — 1, x) (blue circles) and p(N — M+ 1, x) (green squares) in the cosine- 
squared lattice (blue doted line). We clearly see that the probability for particles to 'sit' 
at the minima of the cosine- squared lattice is the highest for N = M atoms (also, the 
probability for the atoms to sit at potential maxima is the lowest for N = M atoms). This 
observation confirms the indication given by SP energy spectrum in Figj2] regarding where 
the pinning takes place in the finite size system. In Fig. [5] (b) we show the same quantities 
for the sine-squared lattice. We see that the signature of pinning is strongest at N = M — 1 
consistent with the single-particle spectrum. 

We see that the energy and the single particle density can distinguish between different 
types of boundary conditions, whereas the GSF is less sensitive to these effects. The GSF has 
an advantage over the energy spectrum in the thermodynamic limit where it dramatically 
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FIG. 4: Single particle density p(N,x) obtained from the many-body ground states of N atoms in 
Vq(x) + Vi(x) potential (a) p(40,x) (blue line) in cosine-squared lattice (red line) and (b) p(40, x) 
(blue line) in sine-squared lattice (red line). For reference we plot the corresponding densities 
po(N,x) of N atoms in the ground state of Vq(x) trap (black doted line). See text for details. 

shows where the pinning takes place for an infmitesimally small lattice amplitude. 
B. Analysis of GSF via 1st order perturbation theory 

In this subsection we study the GSF in the context of the pinning phase transition for the 
box potential via stationary first-order perturbation theory. Unperturbed states are the SP 
basis of Vq(x) i.e ip n (x) = (2/L) sm(nirx/ L). For the moment let us focus on the cosine- 
squared lattice Vi(x) = V\ cos 2 (M7rx/L), which we treat as the small stationary perturbation 
and denote SP basis of Vq(x) + Vi(x) as <p n (x). To first order in the lattice amplitude, the 
single-particle states of the potential V (x) + Vi(x) are 

<f>i(x) oc ipi(x) + a 2M -iip2M-i(x) + a 2M +i'4 , 2M+i(x), (9) 
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FIG. 5: Single particle density, (a) p(40, x) (red cross) and two nearby densities p(39, x) (blue 
circles) and p(41,x) (green squares) in cosine lattice (blue doted line), (b) p(39, x) (red crosses) 
and two nearby densities p(38, x) (blue circles) and /o(40, x) (green triangles) in sine lattice (blue 
doted line). A smaller range of axis is chosen to provide good visibility. See text for details. 

where i G {1, . . . , M — 1, M + 1, . . . , 2M — 1}; the case i = M is treated separately. This 
interval of indices i cover all particle numbers of interest i.e N = 1, . . . , 2M — 1, and the 
criticality region N ~ M is in the center of that interval. The coefficients are given by 
u-zM-i = Vi/[A(E 2M -i ~ Ei)} and a 2M +i = Vi/[A(E 2M+ i - Ei)] where Ei is the SP energy of 
i-th state in the Vq(x) potential. Since Ei = (i/M) 2 En we can write 

M Vi 



(10) 



The coefficient a 2 M+i in Eq. ([9]) can be ignored because of the denominator in Eq. ( ITUj) . i.e., 
for i ^ M 



ipi(x) + a 2M -iip2M-i(x) 



P2M- 



(11) 
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where we have normalized the wave function. We see that perturbation will be most domi- 
nant when i — M — 1 and i — M + 1 . 

For i = M, the first-order perturbation theory gives 4>m{.%) oc 4>m{ x ) + a-iM^zMi.x) , where 
&3M — Vi/[4(E3M ~ Em)] = Vi/32Er is sufficiently small for a weak lattice Vi/Er <C 1, and 
we can write 

Mx)nii> M (x). (12) 

In fact, this relation will hold even for deeper lattices as long as Vi/2>2Er <C 1. 

In order to calculate the ground state fidelity F = | det A| 2 we need to evaluate the matrix 
elements = f ip*(x)(f)j(x)dx [see Eqs. ([3]) and flSJ)]. We first consider the case N < M. 
We use Eqs. ffTTj) and (]12[) to get matrix elements within first order perturbation theory: 

5ij + Cl2M-j5i : 2M-j /10 n 

A a ~ h , | ( 13 ) 

V 1 + \ a 2M-j\ 

where i,j=l,...,N. If N < M, then the second delta term in (fT3|) is zero, and the matrix 



PI) is diagonal A« = (v 7 ! + I^M-i) 2 )" 1 - Thus, the GSF (iV < M) is 

AT Ar 

^ni^-riTTK^F (14) 

Since the coefficients \a,2M-i\ 2 rise quadratically as i approaches M, we understand the 
behavior of GSF when N approaches M from below, which was observed numerically in 
Fig. [1] (a) and (b). In Fig. []Ja) we plot results of Equation f lT^|) (black crosses), agreement 
with exact numerical results is excellent for small amplitudes, i.e, V\ = O.IEr in Fig. []Ja) 
(blue squares), for larger amplitudes agreement is good outside the dip of GSF, while we see 
discrepancies in the dip, resulting from first order perturbation theory are systematically 
lower than exact numerics. 

The case N = M is straightforward due to Amm ~ 1 ; and for N = M the GSF becomes 

Af-i 1 

F ^Ut^i T2-\Amm\\ (15) 

1=1 1 1 

which is identical to the value for N = M — 1 (see Eq. ( fT4l) ). which explains our numerical 
observation. This is an interesting observation. The GSF will decrease when first order 
perturbation is the most effective. We expect it to be the most effective at commensurability, 
N = M. However, the coefficient at N = M — 1 contributes the most in this sense, whereas 
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for iV = M the perturbation on the SP eigenstates (which is reflected onto the many body 
eigenstates via FB mapping) is essentially negligible. 

If we enlarge the box to new size e.g L' = 2L, and we keep the lattice wave vector constant 
k = Mtt/L, we have M' = 2M, and the fidelity dip moves to N = 2M - 1 and N = 2M, 
and decreases in value because the new coefficient clm 1 = 2a m, which gives smaller product 
terms in equation (fl4|) . This explains results of Fig[]Jb) where we vary system size L. 

The first-order perturbation theory also provides an explanation for the influence of the 
phase 0. For example, for the sine-squared lattice, because cos 2 (kx) = 1 — sin 2 (kx), the 
integrals appearing in the perturbation expansion are 

and the coefficients change sign, wave functions differ, but the GSF (TBI) depends on absolute 
squares of these coefficients and is insensitive to this phase. This explains results of Fig. [2] 
for cosine-squared and sine-squared lattice (red circles and blue crosses). For some arbitrary 
phase value between and tt/2, the main difference is that equations (Q and ()12p will 
no longer hold and more coefficients are needed in expansion of (f>i(x) in terms of il>i(x), 
especially for i ~ M, which breaks the symmetry between N = M and N = M — 1 cases. 

Let us finally discuss the N = M + 1 case for the cosine- squared lattice. Matrix A 
acquires two off diagonal elements Am+i,m-i and Am-i,m+i with the following property 

Am+i,m-i — —Am-i,m+i ~ , M+1 , (16) 

a/1 + \a M +i\ 

due to Eq. ([10]) . In this case the determinant of matrix (fT3l) becomes 

M-2 

detA^ j J AhAmm(Am-i,m-iAm+i,m+i — Am+i,m-iAm-i,m+i)- 

i=l 

Due to (TTUD and (USD we have {A M -i,m-iA m +i,m+i - ^m+i,m-i^m-i,m+i) ~ ^ and since 
Amm ~ 1, we finally get that determinant of matrix (|T5|) for M + 1 particles is 

M-2 



det A w ] [ Ah. 



i=l 



From the last relation we see that the ground state fidelities for N = M + 1 and N = M — 2 
particles are approximately equal (see Eq. (JHJ)). In Fig. (Ha) we plot these results for 
N = (M, M+1) particles (black crosses) in addition to GSF for N < M. One could proceed 
to other values of N > M + 1 in the same fashion and analyze the GSF via perturbation 
theory. 
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V. THE PINNING TRANSITION OF THE TONKS-GIRARDEAU GAS IN THE 
HARMONIC OSCILLATOR: GROUND STATE FIDELITY 



In this section we study the pinning transition in the experimentally relevant harmonic 
oscillator (HO) potential [7] 

VM«) = (17) 

We choose parameters following the experiment in Ref. pj], i-e., the atoms inside the trap 
are caesium atoms, 133 Cs. The optical lattice is Vi(x) = Vi sin 2 (kx) with k = 1.887TXo 1 where 
Xo = lfim and the wavelength of the lattice is A = 1064nm. Lattice amplitude and all other 
energies are in units of the recoil energy Er = (hk) 2 /(2m). 

In Fig. |3(a) we plot the ground state fidelity as a function of number of particles for 
different values of the lattice amplitude V\. The harmonic oscillator frequency is u}q/2ti = 
25Hz (similar to the frequency used in experiment u exp /2ir = 22 (3) Hz 7J). We see for all 
amplitudes V\ that GSF first decreases smoothly, in a similar fashion as in the infinitely deep 
box in Fig. d](a) and (b). Then GSF decreases faster until it reaches first minimum at N ~ 50 
after which it develops oscillations with deeper minima. For smaller amplitudes VJ = 0.15Er 
(green circles) and V\ = 0.45-Er (blue x's) the oscillations have a global minimum at N « 60, 
and after that the GSF starts to rise, as expected, but at a slower rate compared to the rate 
of decrease towards the first minimum, in contrast to the square well. For higher amplitude 
V\ = 0.9£'ij(red squares) GSF is effectively zero i.e F ~ for interval of iV's between 
iV ~ 60 and N ~ 70, with the slow increase of average GSF for N above 70. Finally, for still 
higher amplitude V\ = IAEr (black diamonds) (similar to the amplitude used in experiment 
V exp = 1.5(1 )Er [7]) GSF drops to values slightly above zero already for N ~ 45 and after a 
small bump we see that F ~ from iV ~ 60 to N ~ 110; above N ~ 110 GSF slowly rises 
and develops oscillations (not shown) similar to GSF for V\ = 0.9i?^(red squares). 

In Fig. [6](b) we plot the ground state fidelity F as a function of N for different values of 
the harmonic oscillator frequency uq with constant lattice amplitude V\ = 0A5Er. We see, 
as expected, that the pinning transition occurs for larger N, and the fidelity dip lowers. 



We point out that the results of Fig 
with the experimental results of Ref. 



BJa), obtained for 7> 1, are in excellent agreement 
obtained for large but finite 7, where it is said 
that commensurability of superfluid phase and the lattice is best fulfilled when there are 
about iV ~ 60 atoms in the central tube. We see in Fig. E(a) for all amplitudes GSF shows 
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FIG. 6: Ground state fidelity for the pinning transition in a harmonic oscillator potential, (a) GSF 
as a function of number of particles iV for different lattice amplitudes V\ with constant frequency 
ujq/2tt = 25Hz. (b) GSF as a function of N for different ojq with constant Vi = 0A5Er- See text 
for details. 

enhanced sensitivity and strongest decay in the region N ~ 60. 

In order to understand these results, we need to define the commensurability of the 
Tonks-Girardeau gas and the applied optical lattice. This is not straightforward because of 
the inhomogeneous atomic density in the ground state of the harmonic oscillator potential. 
We draw upon the results of Section IIVt where the GSF had a minimal value when the 
unperturbed Nth SP state entering the iV-particle ground state via Eq. (j4j), had the same 
wavelength as the optical lattice. In the case of the harmonic trap (I17j) . the asymptotic 
expansion of SP states i> n { x ) for n ^> 1 is 




(18) 



where ao = \/ h/muo. This provides us with the dominant wavelength of the n-th SP state. 
We estimate that the pinning occurs when 



(19) 



which yields 



N, 



k 2 h 



pinn 



(20) 



2mwo 

for the number of particles where pinning occurs. Eq. ( |20l) is obtained for n > 1; in 



experiments one usually has N > 30. In addition, we stress that Eq. ( 1201) is in agreement 
with Ref. where the pinning transition is said to occur for such iV that the peak density 
of the superffuid phase, obtained with Thomas Fermi approximation for 7 3> 1, is equal 
to the commensurate density n c = 2/ A. For u/2tt = 25Hz our estimate yeilds N } 



pinn 



50, 



which explains the drop in the GSF observed in Fig. HJa). Equation ( 120]) also explains the 
positions of the first minima of GSF in Fig. El^b) since it gives N pinn ~ {38, 52, 75, 130} for 
ujq/2tt = {35, 25, 17.5, 10}Hz, respectively, in fair agreement with exact numerical results. 
Again, as we make the system larger, the minimum value of GSF decreases which is consistent 
with the decrease of GSF at criticality in the thermodynamic limit. 




FIG. 7: Single particle energy spectrum of harmonic trap with the lattice Vq(x) + V[(x), for two 
values of amplitude Vi = IAEr (blue circles) and Vi = 0.45-Er (red triangles). See text for details. 



Besides the GSF, it is instructive also to look into the single particle energy spectrum 
plotted in Fig. [7] for Vq(x) + Vi(x) with the amplitudes Vi = 0A5Er (red triangles) and 
Vi = IAEr (blue circles). We see that for VJ = 0A5Er (red triangles) at n ~ 50 and larger 
values a series of 'gaps' open up in the sense that at some iV values the excitations of the 
TG gas from the ground state cost more energy. 

These results can also be understood simply in terms of commensurability of the SP 
density of the TG gas in the HO potential and the lattice Vi(x). SP density of N particles 
in the ground state of TG gas in Vq(x) potential is po(N, x) = ^2^=0 IVViO^) | 2 - This function 
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is inhomogeneous with maximum value (peak density) at x = (we now ignore the small 
auxiliary oscillations which can lead to the local minimum for x = depending on the parity 
of n = N — 1th state). When we increase N (starting from N = 1), the first part of the 
density po(N,x) that approaches the commensurability condition n c = 2/A, is the central 
part i.e for po(N, 0). Thus the first 'gap', and first minima of GSF, will appear for N pinn with 



property po(N pinn ,0) « n c , in accordance with Ref. jfi] and Equation ( 120]) . Adding more 
particles, i.e iV > N pinn , leads to po(N, 0) > n c but now in some regions left and right from 
x = 0, the density becomes commensurate with lattice i.e po(N, —d) ~ n c and po{N, d) m n c , 
for some d > 0, and pinning still occurs, i.e, additional 'gaps' in SP spectrum are present 
and GSF still lowers. The distance d increases with the increase of N and commensurability 
condition is satisfied for regions further towards the edges of the trap Vq(x). The fraction of 
the atomic cloud commensurate with the lattice gets smaller and as a consequence the GSF 
slowly increases. Oscillations are present due to many fine details such as the interplay of 
symmetry of the lattice and the symmetry of the i/jn(x) states. 

VI. THE LOSCHMIDT ECHO AND OUT OF EQUILIBRIUM DYNAMICS: SUD- 
DEN QUENCH WITH OPTICAL LATTICE 

In this section we look for signatures of the pinning transition in the non-equilibrium 
dynamics of TG gas, initially in ground state of some trapping potential Vq(x), after optical 
lattice Vi(x) is suddenly turned on. This is an example of a sudden 'quench'. Before the 
quench with the lattice potential, the gas is in the equilibrium ground state |\&o) of Hamil- 
tonian H Q , see ([T]). At t = we suddenly turn on the lattice potential Vi(x) and an out- 
of-equilibrium many body state |$(t)) = exp[— i(Hq + Vi{x))t / fi\\^ q) starts to evolve, where 
|$(0)) = |^o) is the initial condition. We would like to develop a quantitative understanding 
of the out of equilibrium dynamics which is encoded in this state. Conveniently, the mapping 
theorem also holds for time dependent states, and the quenched state \$(t)) can be con- 
structed using a Slater determinant of time evolving single particle states such that |3>(t)) = 
-t=j T[i<i<j<N s & n ( x i ~ x j) det^ =1 [ijj n (xj,t)]. The single particle states ip n (xj,t) are out of 
equilibrium, and obtained by solving ihd t ip n (x,t) = \—h 2 /{2m)d1 + Vq(x) + Vi(x)]i/} n (x,t) 
with initial conditions ip n (x,0) = ip n {x) where ip n {x) are the initial single particle states 
(SP) which are used to construct the unperturbed ground state |^ ) i- e SP states of V (x) 
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potential. 

A prototypical quantityto calculate for a system perturbed out of equilibrium is the 
so-called Loschmidt echo 15j, which is defined as 



L(t) = \(^ \eMiH t/h)ex V [-i(H + Vi(x))t/h}\^ }\ 2 . 

It is a measure of the sensitivity of the system to the quench protocol, which in this case is 
simply the application of the external lattice potential to the initial Tonks gas equilibrium 
state. Despite it's mathematical simplicity, it conveys a great deal of information about the 
many-body system under scrutiny, such as universal behavior at criticality jg] and important 
information on the thermalization of observables. Closed formulas for the echo are, in 
general, very difficult to obtain. For a Tonks-Girardeau gas, the Loschmidt echo was recently 
computed in a relatively straightforward way 17 1. 

Since Hq\^q) = fio|^o)j where Oq is ground state energy of TG gas in Vq{x) trap, we get 



Lit) = |(vl> |expH(if + ^(x))t/n]|vl> )| 2 = \(^ \<$>{t))\ 2 . (21) 

Relation (121]) shows that in our case the Loschmidt echo (LE) is equivalent to the survival 
probability i.e probability that system will be in the initial state at the time t after the 
quench. We will interchangeably use the terms Loschmidt echo and survival probability. 

The Fermi-Bose mapping theorem is valid for time dependent wave functions thus LE 
can be written in a form convenient for calculation, analogous to the calculation of the static 
fidelity, 

N N 

L{t) = i M /&r- ^ nc W (^«, o) n ^w^'*)! 2 

CTl i = l CT2 J=l 

= I det A(t)| 2 (22) 



and A(t) is the time dependent matrix containing overlaps between static SP states of the 
Vq(x) potential i.e ipi(x, 0) and SP states if)i(x, t) evolved in perturbed potential Vq(x) + Vi(x) 



Aij(t) = / ^*(^^)^ j (x,t)dx. 



(23) 



Equations ( 1221) and ( 1231) were recently used to study the long time behavior of many-particle 
quantum decay 20|. The LE of the TG gas is formally equivalent to the corresponding echo 
for a gas of non-interacting fermions li| . The Loschmidt echo of one-dimensional interacting 
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3ose gases was recently related [17| to series of experiments [2lM23| and theoretical studies 
24j-|29| on interference between split parallel ID Bose systems. 



A. Infinitely deep well 

In this subsection we use (1221) and fT2"3"l) to explore the LE of a TG gas after a sudden 
quench with optical lattice Vi(x), i.e at t = we suddenly turn on the lattice and leave it 
on. Before the quench TG gas is in ground state |\l/o) of infinitely deep well flH]) potential 
Vo(x). In this subsection we use same units and parameters as in section HVl 

In Fig. [S] we show Loschmidt echo L(t), as a function of time for different num- 
bers of particles N obtained with exact numerical evolution, for optical lattice V\{x) = 
0.55Eji cos 2 (47rxa; ( ^ 1 /3). Size of the well L = 30xq. Number of lattice wells is M = 40. 
That set of parameters is the same as the one we used for ground state fidelity F in Fig. 
[D[a) denoted with red circles. We see that properties of F are reflected in LE. Values of 
LE goes in pairs, i.e, curves are approximately equal for N = M — j and N = M + j — 1 
particles, with j — 1, . . . , M — 1. Decay of Loschmidt echo is strongest and fastest close to 
criticality, i.e, for N = 39 and N = 40 particles (red circles and red solid line). In addition, 
the oscillations get slower as we approach criticality. 
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FIG. 8: Loschmidt echo L(t) for different particle numbers N in the infinitely deep box potential. 
Optical lattice Vi{x) has 40 minima, i.e, M = 40 lattice wells. LE reflects properties of GSF, i.e, 
decay is strongest for N = (M — 1,M), and Loschmidt echo appears in pairs. The frequency of 
revivals decreases as we approach criticality. See text for details. 

In order to understand the results of Fig. [8] we use relations (|T0|) - (fT2|) to write expansion 
of the unperturbed SP states of Vq(x) in terms of the perturbed SP states of V (x) + Vi(x) 
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&{ x , 0) « Ux) = ^Mz^MzM (24) 
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2 exp(-iE l 2M _ i t/H) 
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For i = 1, . . . , M - 1, M + 1, . . . , 2M - 1. In case of i = M we get 

1>m(x,0) ™ Mx)- (25) 
Using relation fl2^|) we get time evolution of after the quench 

a/1 + |a 2 M-il 

where are SP energy's in Vo(x) + VJ(x) potential. Now we use (p4l) - (|26l) in Equations (|23|) 
and (122|) in order to obtain insight in to the behavior of the Loschmidt echo. The following 
analysis is very similar to analysis done in subsection IIVBI for the ground state fidelity. 

First, we consider the case when the number of particles is less than the number of wells, 
i.e N < M. In this case matrix fl23l) is diagonal with 

and since L(t) = | det A(t)| 2 = det [A(t)A*(t)] we get 

N 

Ht)*l[\Mt)\ 2 , (27) 

1=1 

where 

lA , A g „ 1 + 2|a 2Af - t | 2 cos [(E l 2M _ { - Ept/H] + \a 2 M-i\ 4 

lAiWI — — r~2 • (/°J 

(1 + \a 2 M-i\ 2 ) 

If we use expression Q14p for the ground state fidelity (F) we get for the Loscmidt echo of 
N < M particles 

N 

L(t) « F 2 J](l + 2|a 2A ,_ l | 2 cos [{E l 2M ^ - E^t/h] + |a 2M -i| 4 ), (29) 
i=i 

we see that ground state fidelity is incorporated in LE by construction i.e., that is why the 
LE reflects it's properties. Since coefficients \a 2 M~i\ 2 grow quadratically as i approaches M 
the most dominant cosine term in ( 129]) is for i = N which leads to the conclusion that the 
dominant frequency of revivals ur for LE of N < M particles is simply related to the SP 
energy of Vq{x) + Vi(x) potential through the relation 

uJr{N) = . (30) 
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Now consider the case N = M . Going from N = M — 1 to N = M particles, the matrix 
fl23|) remains diagonal due to relation f )25|) . we simply add AMuif) ~ exp(—iE l M t/h) on the 
main diagonal, and since |^4mm(^)| 2 ~ 1, we get that the LE for N = M is the same as the 
LE for N = M — 1 particles, in accordance with Fig. [5] (red solid line for N = 40 and red 
circles for N = 39 particles). 

Now we proceed to the N = M + 1 case. In this case the matrix (1231) obtains two off 
diagonal elements 

a (+ \ A ,+\ „ -aM+iexp(-iE l M+1 t/h) - a M -i exp(-zE^ 1 t/^) 

Am+\,M-1\J>) — AM-l,M+ll£J r T~ i 19 ) l^-U 

1 + \a M +i\ 

where we used (1241) and |om-i| = |a A/+1 1 , see fTTUj) . The determinant of the matrix f l23|) 
has two terms, the product of diagonal elements and a term arising from two off diagonal 
elements 

M-2 

det A(t) w JJ A^AmmO) [^M-i,JVf-i(t)AM+i,M+i(t) - A M+ i,M-i(i)^-i,M+i(t)] • 
i=i 

It can be shown that 

[Am -i,M-i{t) A M+ltM+ i(t) — A M+ltM _i(t)A M _i tM+ i 
which together with A MM {t) ~ exp(—iE l M t/h), yields 

M-2 

L(t) = |detA(t)| 2 ^ J] |^(t)| 2 . 



8=1 



We conclude that the Loschmidt echoes for N = M+l and N = M—2 are approximately the 
same. We see that the same pattern emerges as with ground state fidelity. The Loschmidt 
echo values come in pairs, i.e, it is approximately the same for N = M—j and N = M+j — 1 
particles, where j = 1, . . . , M — 1. Due to this pattern, we can use Eq. (1301) to get the 
dominant revival frequency of Loschmidt echo for other particle numbers iV > M, i.e we 
can write 

u R (M + j-l)Ku R (M-j), (32) 

where j = 1, . . . , M — 1 and ujr(M — j) on the right side of Eq. ( |32i) is given by Eq. ( 1301) . In 
order to check the quality of these relations, we plot in Fig. [9] the dominant revival frequency 
u)r(N) obtained via (|30|) and (1321) for N = 1, . . . , 2M — 2 particles (for parameters used here 
M = 40), together with the most dominant frequency ufft{N) obtained with the Fourier 
transform of LE; the agreement is excellent. Parameters used to calculate the LE are the 
same as in Fig. El 
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FIG. 9: Dominant revival frequency of Loschmidt echo obtained with Eqs. (|30p and (|32p (blue 
circles) and with the Fourier transform of Loschmidt echo ujfft (red crosses). Parameters used 
for calculation of LE are the same as in Fig. See text for details. 

B. Harmonic oscillator 

In this subsection we use Eqs. (1221) and (1231) to explore the LE of a TG gas following 
a sudden quench with optical lattice Vi(x); before the quench the TG gas is in the ground 
state |\l/o) of harmonic oscillator potential Vq(x) (fTTj) . In this subsection we use same units 
and parameters as in Section PVT. 

In Fig. [10] we show the Loschmidt echo following the quench with optical lattice Vi(x) = 
0.45-Etj sin 2 (1.887ncQ x); the system was initially in a ground state of the harmonic oscillator 
potential with cuq/2-k = 25Hz. Fig. [TU](a) is for iV < N pinn and Fig. [Hjb) is for N > N pinn , 
where N pinn is given by Eq. (12~DT) (for parameters used here N pinn pa 52). We see that 
properties of GSF (see Fig. [6](a) blue crosses) are reflected in the LE; the average values of 
the LE are lower for lower GSF. This is a general observation. However, the details of LE 
dynamics (such as the dominant revival frequency) depend on trapping potential. 

In Fig. [TTJ we illustrate the dominant frequency of the LE (from Fig. ITOT) obtained with 
the Fourier transform (ojfft blue crosses). We see that ujfft is constant for iV < N pinn , and 
it starts to behave irregularly for N > N pinn . We have found numerically that the regular 
behavior for N < N pinn occurs because ip n (x, 0) can be well approximated with 

ip n (x, 0) ~ A(j) n _ 2 (x) + B(f) n (x) + C(p n+2 (x), 

where coefficient B is always the largest in magnitude; this yields 

ur(N < N pinn ) ^ 2uj (33) 
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FIG. 10: The Loschmidt echo L(t) for different particle numbers N. a) LE for particle numbers 
up to N < N p i nn . b) LE for N > N p i nn . See text for details. 

which is also plotted in Fig. [TTJ For N > N P i nn there are many coefficients in the expansion 
of ip n (x, 0) in terms of (j) n (x) which contribute on equal footing and simple relation (|33|) does 
not hold anymore. 
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FIG. 11: Dominant revival frequency of Loschmidt echo obtained with Eq. (|33[) ur (red circles) 
and with the Fourier transform of the Loschmidt echo oofft (blue crosses). The parameters used 
to calculate the LE are the same as in Fig. [Tol See text for details. 
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VII. CONCLUSION 



We have studied the pinning quantum phase transition in a Tonks- Girardeau gas, both in 
equilibrium and out-of-equilibrium, using the ground state fidelity and the Loschmidt echo 
as diagnostic tools. We have found, both numerically and analytically (within first order 
perturbation theory), that the ground state fidelity can individuate the region of criticality. 
The ground state fidelity defined in Eq. (jSJ) has a dramatic decrease when the atomic density 
approaches the commensurate density of one particle per lattice well. This decrease is a 
signature of the pinning transition from the Tonks to the Mott insulating phase. We have 
found that the ground state fidelity of the TG gas in an infinitely deep well potential can 
be insensitive to finite size effects. The GSF for N = M — 1 and N = M particles (M 
denotes number of lattice wells) in cosine-squared (sine-squared) lattice are the same, while 
the single particle energy spectrum and density show that pinning actually happens for 
N = M (N = M — 1, respectively). The GSF has an advantage over the density and single 
particle energy spectrum in the thermodynamic limit, where it dramatically shows where 
the pinning takes place for an infinitesimally small lattice amplitude. We have studied the 
applicability of the fidelity for diagnosing the pinning transition in experimentally realistic 
scenarios. Our results are in excellent agreement with recent experimental work 7]. 

We have found that the GSF in a harmonic oscillator potentials shows enhanced sensi- 
tivity in a broad region of particle numbers N > N pinn (where N pinn is defined in Eq. (1201) ); 
at N ~ Npi nn GSF has a faster decay and for larger N develops oscillations. This behavior 
is related to series of 'gaps' opening at n ~ N pinn in the single particle energy spectrum of 
the total potential (harmonic oscillator plus optical lattice). 

In addition, we have explored the out of equilibrium dynamics of the gas following a 
sudden quench with a lattice potential potential. We have showed that all properties of the 
ground state fidelity are reflected in the Loschmidt echo dynamics i.e., in the non equilibrium 
dynamics of the Tonks-Girardeau gas initiated by sudden quench of the lattice potential. 
The average value of the Loschmidt echo is lower for lower values of ground state fidelity, 
regardless of the details of the trapping potential. Details of the Loschmidt echo dynamics, 
such as dominant revival frequency, depends on the type of trapping potential. We find 
regular behavior of revivals for all relevant particle numbers in infinitely deep well potential 
i.e, frequency's get lower as we approach criticality and can be calculated simply from single 
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particle energy spectrum of the total potential (infinitely deep well plus optical lattice). In 
the harmonic oscillator potential, the dominant frequency of revivals behaves in a regular 
way. It is a constant, approximately equal to 2uo (uq is frequency of harmonic trap), until 
a series of 'gaps' open in the single particle energy spectrum of the total potential. 
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